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ABSTRACT 


We calculate the distribution function of astronomical objects (like galaxies 
and/or smooth halos of different kinds) gravitational hel ds due to their t i dal in - 


teraction. For that we apply the statistical method of 


Chandrasekhar! (119431) . 


used there to calculate famous Holtzmark distribution. We show that in our ap¬ 
proach the distribution function is never Gaussian, its form being dictated by the 
potential of interaction between objects. This calculation permits us to perform 
a theoretical analysis of the relation between angular moment um and mass (rich¬ 


ness) of the 


galaxy clusters. To do so , we follow the idea of 


Catelan fc Theuns 


(1l996l) and iHeavens fc Peacock! (119881) . The main difference is that here we re¬ 


duce the problem to discrete many-body case, where all physical properties of the 
system are determined by the interaction potential V{rij). The essence of reduc¬ 
tion is that we use the multipole (up to quadrupole here) expansion of Newtonian 
potential so that all hydrodynamic, ’’extended” characteristics of an object like 
its density mass are ’’integrated out” giving its ”point-like” characteristics like 
mass and quadrupole moment. In that sense we make no difference between 
galaxies and smooth components like halos. We compare our theoretical results 
with observational data. 


Subject headings: galaxies: clusters: general 
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Introduction 


The problem of galaxies and their structures formation is one of the objectives of 


formation 

Peebles 

1969 

: Sunvaew & Zeldovich 

1972; 

Zeldovich 

1970 

Shandarin 

1974; 

Dekel 

1985; 

Wesson 

1982; 

Efstathiou & Silk 

1983). 


Doroshkevich 


1973 


important. This is because the new scenarios are essentially the modihcations of the 
old ones and can be classihed according to the classical ones . Revised and improv ed 


structure formation scenarios can b e found in va r ious papers flLee fc Pen 


2000 


2001 


2002 


Navarro et ah ^ 

>004; 

VIo et al. / 

>005 


Bower 

2005; 

Trnii 

lio et al. 

2006; 

Brook e 

al. 

2008 


Vera-Ciro et al 

2011 

: Paz et al 

2008 

Shandarin et al.l 

2012; 

Codis et al. 

2012; 

Varela et al. 


2012 


Giahi-Saravani fc Schafer 


20141) . The crucial goal is to discriminate among different 


models of galaxies and their structures formation. The main controversy here is how 
galaxies acquire their angular moments, which yield subsequently the moments of galaxy 
clusters. 

Presently the commonly accepted is spatially flat, homogeneous and isotropic ACDM 


model of the Universe. In such model, the structures were formed 


rom 


he primordial 


197ol: Sunvaew & Zeldovich 

1970 

). The most popular galaxy 

brmation scenario, so-called 

hierarchic clustering model 

^Doroshkevich 1970 

Dekel 

19§5; 

Peebles 

1969 

) is 

)ased on 

this assumption. The numerica’ 

simulations (I 

Bond et al. 

1996; 

Sprineel et al. 

2005: 


van de Wevgaert fc Bond 


2nn8a 


bl) conhrm that such mechanism could be realized in the 


Universe. In this mechanism, the large scale structure can appear from bottom to up as a 
consequence of gravitational interactions between galaxies. This means that galaxies are 
formed at the beginning with subsequent merger into larger clusters (structures). In this 
case, the galaxies spin angular momenta arise as a result of interaction with their neighbors. 
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Original version of this model claims that the initial orientation of galaxy spins should 
be random. However, it had been shown later, that hierarchic clustering model admits 
naturally so-called Tidal Torque mechanism. In this mechanism, galaxies have their angular 
momenta aligned due to the coupling between the protogalaxy region and surrounding 


structure . The galaxy rotation in 


galaxies (jWesson 


1982 


White 


his mechanism is due to tidal interaction between 
le ideas of 


1984^ based on t 


Tidal Torque scenario is presented by 


Hovld (jl95lf) . The review of 


Schaeferl (120091) . Note that while originally the Tidal 


Torque mechanism considers completely random distribution of spin angular momenta of 


galaxies, it has been shown later t. 

lat the local tidal shear tensor can cause a lo 

cal alignment 

of their rotational axes ( 

dubinski 

1992; 

Catelan & Theuns 

1996 


Lee & Pen 

2000, 

2001 


2002; 

Navarro et al. 

2004 

). On the other hand, some authors like 

Brook et al 

(200 

8), still 


argue that there is misalignment of angular momenta in the hierarchical clustering model. 

Different scenarios make different predictions concerning distribution of their angular 
momenta and esp ecially about orientation of galaxies in structures (for review see 


Godlowskil (l2011al) ). Of course the hnal test of a given scenario is to compare its predictions 


with observations. It is possible to conclude that the observed variations in angular 


(Romanowskv & Fall 

2012; 

Joachimi et al. 

2015) 


The investigations of galaxies alignment show that it depends on the mass of the parent 
structure. Generally, the groups of galaxies and the small galaxy clusters reveal almost 


no alignment, while we observe such alignment for ri ch galaxy clusters anc 


The alignment increases with mass of t 


Godlowski et al. 


2010 


Godlowski 


2011a 


le structure f Godlowski et al. 


2005 


superc 


201211 . 


Arval 


usters. 


2007 : 


Note that it is commonly agreed that there is no evidence for rotation of the groups 
and clusters of galaxies. That implies that such structures do not rotate (for example 
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Regos fc Gelle: 


Hwang fc Led (120071) ). see however 


( 1989 ): Diaferio fc Geller ( 1997 ): Diaferio ( 1999 ): Rines et al. ( 2003 ): 


Kalinknv et a 


this context, especially important is the result of 


(120051) for an opp osite opinion. In 


Hwang fc heel (120071) who examined the 


dispersions and velocity gradients in 899 Abell clusters and found possible evidence for 
rotation in only six of them. Thus any non-zero angular momentum in groups and clusters 
of galaxies should arise only from possible alignment of galaxy spins and stronger alignment 
mean larger angular momentum of such structures. 

The aim of the present work is the theoretical analysis of the influence of tidal 
interaction between objects like galaxies, their clusters as well as smooth halos on their 
gravit ational held distribution. For that we use the statistical method of 


Chandrasekhar 


(119431) . To apply our result to observable quantities, we calculate th e distribut i on fun ction 


of the angular momenta. We do that for linear (corresponding to 


Zeldovichl (119701) 


approximation in displacement held) and nonlinear regimes of huctuation growth. As in 
our method, the parameters of galaxy ensemble like their masses, radii and volumes enter 
problem as parameters, our calculation permits to trace possible relation between angular 
momentum and mass (richness) of the galaxy clusters. The above statistical method 
reveals the fact, that in the stellar systems, the derived distribution function cannot be 
Gau ssian but rather be l ongs to the family of so-called ’’heavy-tailed distributions”, see. 


e.g., 


Kapur fc KesavanI (119921) for details. Moreover, choosing the cosmology on the base of 


corresponding Friedmann equation, our result permits to trace the time evolution of the 
distribution function of angular momenta and its mean value L. In our approach, we can 
also derive the well-known empirical relation between mean galaxy cluster moment L and 
its mass M, L ~ . 


The paper is organized as follows. In Sections 2 and 3 we, based solely on the 
quadrupole (tidal) interaction potential between astronomical objects, calculate our 
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universal (i.e. independent on the details of Lanrangian or Eulerian spaces) distribution 
function of gravitational fields f{E). In the Section 4, on the base of the function f{E), 
we calculate the distribution function of angular momenta f{L) both in the linear and 
nonlinear Lagrangian approximation. We emphasize that as f{E) does not depend on 
the Eulerian or Lagrangian picture, it can be used to calculate the distribution of any 
quantity (like momentum) of the astronomical objects (not only galaxy clusters but smooth 
component like halos as well) in any (linear or nonlinear) regime of fluctuation growth. 


2. General formalism 


Unfortunately, the spin angular momentum is known only for very few galaxies and 
structures. Therefore, instead of the angular momentum by itself, the orientation of galaxies 
in each cluster is usually studied. This is the reason that here we are interested primarily 
in the absolute value of galaxies angular momenta. 


We represent a matter (both luminous and dark) as the Newtonian self-gravitating 
fluid embedded in the Universe obeying corresponding Friedmann equation. To obtain the 
tidal (i.e. shap e distorting) interaction between the astronomical objects, we, similar to 


fjPoisson 


19981 ). perform the multipole expansion of the Newtonian interaction potential 


between flu id elements. T runcating this expansion on the quadrupole terms, in the spirit of 


the article fjPoisson 


19981) . we write the Hamiltonian (total classical energy) of interaction 


between galaxies in the form 


H = -G ^ QiiTijV(vij) = -GM"^ '^pimiTrijV (r^), 
ij ij 

1 3 cos^ 9ij — 1 


( 1 ) 


Vir,,) = 




where G is the gravitational constant, M = is the total mass of the ensemble. 

Pi = Qi/Qi is the quadrupole moment of the Uth galaxy, ruj is its mass, rij = |rjjj. 
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Tij = Tj— Ti is a relative separa tion between c entres of galaxies and N is their number. The 


expression for Qi has the form fjPoisson 


199811 


Qi= pi(x)|xpP2(s • x)ci^a;, P2{x) = - 1), 

J Vi 

where 14 is a volume of i-th galaxy, pj(x) is a density of its mass. 


( 2 ) 


Note, that expression ([T]) generalizes the two-particle result of Ref. fjPoissonI 1199811 on 
the ensemble of N objects, splitting the interactions between them in pairwise manner. 
Such splitting is customary in condensed matter physics where the interacting many-body 


ensemble is represented by the sum of all 


(for example 123=12-1-13-1-23), see Ref. f Mailis 


possib 


e pai r interactions between particles i and j 


200011 and references 


of magnetic systems. The same procedure is also used in astronomy (jPeacock 


herein in the context 


19991). 


The physics of the interaction ([T]) is the following. Under the influence of the interaction 
with other astronomical objects, the shape of a given i-th object changes, which inflicts the 
variation of its density held Pii^x.), Eq. ([2]). As the objects are situated randomly and have 
random shapes, the mass rrii and quadrupole moment Qi of an object (like galaxy, galaxy 
cluster or smooth component) will vary randomly. This generates the random variations of 
the gravity held Equad from these quadrupoles. Latter held is a gradient of the potential ([I]) 
(divided by mass m) and has the form 


E,w(r) = \ = ^ 

where is the unit vector in the direction of the radius-vector. 


(3) 


As we cannot solve the random many-body problem ([T]) exactly, diherent approaches 
had been used for its approximate solution. The simplest approach in condensed matter 
physics is so-called mean held approximation, where the huctuating (electric, magnetic, 
elastic in condensed matter) held, acting on the spec ihc particle is substituted by some 


ensemble averaged, mean held, see, e.g. flMailis 


2000l) . This approach does not take into 






















account the particle clustering (i.e. two-, three-, four-particle cluster in the held of the 
rest of likewise clusters) and is not suitable to describe the systems with disorder, which is 
also our case. In the disordered cas e the adequate appro ach is statistical method similar 


to that applied by Chandrasekhar (jChandrasekhar 


19431) to describe the huctuations of 


the force acting on a specihc star from the rest of stellar system. In this method, the 
distribution function of random gravitational held is introduced so that any observable 
quantity of an ensemble (like orbital moment, energy etc) can be expressed by averaging 


the corresponc 


flStephanovich 


ing s in gle-particle quantity with the a b ove distribution function, see Re fs. 


19971) . flSemenov fc Stephanovich 


2002 '). ( Semenov fc Stephanovich 2003 ) and 


references therein, where the statistical method had been applied to disordered solids. 


3. Distribution function of quadrupolar fields 


According to sta tistical method, very similar to that from Chandrasekhar 


Chandrasekhar 

L943 


(see 

Stephanovich 

(1997) 


19431). the distribution function o 


rando m quadrupolar helds reads 


/(E) = <5(E-E,), 


(4) 


where 3(x) is Dirac 6 - function, Ej is given by the expression (|3]) and bar means the 
averaging over spatial (and any other possible) disorder. We note that the above statistical 
method ha s frequently been us ed to describe the physical properti es of disordered solids like 


29931) . 


dielec trics flStephanovich 


19971) and/or magnetic semiconductors f|Semenov fc Stephanovich 


The physical meaning of distribution function (j4]) is as follows. If we do not have any 
randomness in the galaxies ensemble (for instance, all galaxies are similar) the distribution 
function is just delta-peak, centered in the corresponding held Ej. If we have a disorder in 











































a system, the averaging broadens this delta-peak, giving rise to ” bell-shaped” continuous 
probability distribution. During this averaging the index i will be ’’swallowed” as we 
average over the stellar ensemble. 

To perform the prescribed averagings, we pass to the integral representation of Dirac 6 
- function to obtain 

_ 1 _ 

/(E) = <5(E - E,) = / e^^Pe-^P^^dp^dpydp,. (5) 

J —oo 


The explicit averaging in eq. (E]) is based on the fact that the mass and quadrupole 
moment of the galaxies in the volume V obey the uniform distribution with probability 
density equal to 1/D. This is equivale nt to the fact that th e fluctuati ons of above gala xies 


parameters obey Poisson distribution fj Chandrasekhar 


I 943 I) . see also ( Zolotarev Ill9861) for 


purely mathematical treatment of this question. In this case the averaging for single galaxy 
yields 


[g-ipEU ^ 


47rD 


( 6 ) 


'V 


while for the N galaxies ensemble the correspond ing average equals to e = [e 


i.e. the single-galaxy average, raised to power N fj Chandrasekhar 


1943h . 


We note here that the above averaging procedure im plies that the astronomica 


are similar to each other. However, it can be shown (see 


Semenov fc Stephanovich 


objec ts 


and references therein), that such approach takes into account the pair clusters of galaxies 
or other spatially disordered constituents in the systems other then stellar. The higher 


order c 


(jZiman 


listers (three, four etc objects) can be taken into account along the lines of Ref. 


19791) which would require to solve the chain of kinetic equations for n = 3,4 


etc bodies distribution functions. On the other hand, since the above procedure accounts 
exactly for pair clusters, the many-body clusters can be considered by the splitting them 
into corresponding pairs. This means that our approach takes also the many-particle 
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clustering into account. However, in the future, we are going to elaborate the above n - 
bodies procedure and compare it to our present results. 


We have hnally 


e = lim 

N^oo,V—>-oo 


1 N 


AtiV 


g-ipE(r)^3^ 


'V 


= lim 

N^OO,V^OO 




U{r)d^r 


’V 


N 


= lim <(1-/ 1 — 17(r) 

Af—>-00,V— [ 1/ Jy 


N 


(Tr ^ = lim < 1 — 

V—>oo 


n 

N 


'V 


1 - U(v) 


N 

d-^r^ = 


= exp 


—n / [l — H(r)](i^r 


'V 


C(r) ^ i- 

dvr J pE{r) 

where E{r) is given by Eq. ([3]). 


(7) 


The expression for U{r) signihes the averaging over random angle between vectors p 
and E. Also, we denote in ((Tj) A^ to be number of galaxies so that n = N/V = const is a 
density of galaxies in a given volume. 

Combining the equations (ISl) and ([7]), we obtain following expression for the distribution 
function 


/(E) = 


giEp-F(p)^3^^ F(p) = n 


1 - 


sin pE{r) 
pE(v) 


d-^r. 


( 8 ) 


(27r)3 jv 

We see that E[p) is indeed the characteristic function for random gravitational helds 
distribution. Note that the galaxy clustering can be better considered if we assume that 
the density of galaxies is not a constant but rather n = rafr). Then, similar to the ca se of 
disordered magnetic semiconductors [see Eq. (6b) of Semenov fc Stephanovich 12003 1]. the 
function ?7,(r) should be put under the integral sign in ([8]) to give 

sin pE(r 


F(p) = / n(r) 
Jv 


1 - 


pEir) 


dr. 


( 9 ) 


If we specify the empirical dependence n(r), the characterstic function E[p) can be 
calculated only numerically. 
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Below we shall calculate this function without any assumptions analytically for the 
simpler case n = const. One more way to consider the effects of clustering is to account for 
inhomogeneous distribution of masse s (and/or quadrupola r moments) in the ensemble. This 


can be done along the lines of Ref. fj Chandrasekhar 


194311 ■ where the distribution function 


of masses r(m) had been introduced. In our case, however, the situation will be not that 
simple as we are dealing with more complicated object like quadrupole moment of an object. 
In the present publication we do not consider this effect especially in view that there is 
large ambiguity in determination of function r(m) from the astronomical observation data. 
However, we are going to incorporate the dependence r(m) in our consideration in future. 

We note here that distribution function /(E) ([8]) is by no means Gaussian. We 
will show that function (|8]) does not admit Gaussian limit so that for problem of 
interaction between gravitational quadrupoles (and actually higher order terms like 
octupoles etc) the distribution function of gravitational helds and angular momenta is 


never Gaussian. To demonstrate 
matter systems fferroelec t rics in 


Tat, we note t 


Stephanovich 


rat it had been shown earlier for condensed 


19971 ) and magnetic semiconductors in 


Semenov fc Stephanovichl (120031) ) that Gaussian limit corresponds to large density n —oo 


of electric dipoles (ferroelectrics) or spins (magnetic systems). 


It had also been shown by 


Stephanovichl fjl997l ) and 


Semenov fc Stephanovichl (120031 ) 


that limit n —)■ oo corresponds to small Fourier variable p in the equation ([8]). This means 
that to obtain the Gaussian limit of the distribution function ([8]), we should expand its 
characteristic function F{p) in small p. In hrst nonvanishing (Gaussian) approximation in 
small p this procedure yields 

.,2 r 


Fg(p) = 


np 


E\r)Fr 


( 10 ) 


'V 


To derive equation dTOj) . we use the expansion sinx/x ~ 1 — + ..., valid at x —0. The 
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explicit rewriting of ffTOj) with respect to ([3]) yields 


1 /-OO fjj. 

-^g(p) = / (3 cos^ P — 1)^ sin PdP / —, 

3 Jo Jo ^ 


( 11 ) 


where the last integral over r is divergent at small r. We note here that in the solids 
the Gaussian limit exists (i.e. the integral (ITT]) becomes convergent) due to presence of 
short range terms ~ exp(—r/r^) (r^ is so-called correlation radius, defining the range 
of interaction) in the interaction potential between dipoles or spins. As the exchange 
interaction between spins is always of short range, this feature is peculiar to magnetic 


systems, see ISemenov fc Stephanovich] (120031) and references therein. This means that due 
to long-range quadrupole-quadrupole interaction between galaxies, ~ r“^, the distribution 
function of their helds and/or orbital moments (see below) is never Gaussian. 


The explicit expression for F{p) reads 


r-27r 


F{p) =n r'^dr / sin 9d9 / dp 


sin pEo(3cos2>9-l) 
__ 

p£^o(3 cos^ 6 — 1 ) 


Jo Jo Jo 

We hrst perform the integration over r in (IT^ . For that we denote h/r^ 
{h = pEo{3 cos^ P — 1)) to obtain 


= X 




Q ) 4'“' ./n 

The last integral can be calculated analytically to give 

dx (^ sinx 


dx /^ sm X 


X 


'0 


X 


7/4 


1 - 


X 


2^,2-V2r(-) ^ 1.05711302, 


where r(a;) is T - function (jAbramowitz Sz Stegun 


1964 ). This finally yields 


K = 0.264278|5|^/^ = A\b\^/^. 


( 12 ) 


(13) 


(14) 


With respect to (fT4|) we have 


F{p) = 2TmA{pEQ)^^^ / sin6*(i6'|3cos^6* — Ip/^. 

Jo 


( 15 ) 
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The auxiliary integral Ki can be calculated numerically to give 

Ki = j dx\3x‘^ - = 2 ^ dx{3x‘^ - 1 )^/^ = 1 . 581940777 . 

With respect to ffTOj) . F{p) assumes the form 

F(p) = 2iTn ■ 0.41807255 ■ (pE„f'\ 


(16) 


(17) 


Substitution of ffT7|) into ([HD gives 


'OO 

p sin pE dp, 


2irn ■ 0.41807255 • . 


(18) 


a 


The expression 0181) is the main theoretical result of the present paper, constituting the hnal 
answer for distribution function of gravitational helds moduli. It is seen that distribution 
function (1T8|) depends parametrically on the galaxies density n, as well as on average galaxy 
quadrupole moment. To the best of onr knowledge, neither distribution function (1T8|) nor 
its explicit dependence on n and Q has been known previously. The integral OTSD can be 
calculated only numerically. 


The normalization condition for distribntion function fll8p looks like 


dvr/ E‘^f{E)dE=l. 

Jo 


(19) 


We check explicitly 

1 = 471 I E^f{E)dE = 


OO r\ 

2 ^ 


OO poo 

-I E sin pEdE pdpe-^p"^\ 
TT Jo Jo 


With respect to the relation EsinpEdE = —7i5'{p), the integral / can be calculated 


poo J 

I = —2 p 6'{p)dp = — p e 

Jo dp L 


-«p3/4 


= 1 . 


Jp=0 


( 20 ) 
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To derive the result 
function 


, we use the following identity for n -th derivative of Dirac 5 


f {x)5^'^\x)dx = — 


^5^'^-^\x)dx. 

dx 


We dually mention the diff erence of our approach t o the problem of angular momen ts 


(Gatelan & Thenns 

1996 

) and ( 

Heavens & Peacock 



The main difference is that we consider the discrete many-body problem, stemming from 
multipole expansion (up to quadrupole here) of Newtonian interaction potential between 
duid elements. In such approach, all physical properties of the system are determined by 
the interaction potential V{rij) ([T]). In the case of disorder, the distribution function of 
random gravitational delds is also completely determined by the form of potential V{rij). 
To derive the distribution f unctio n of random gravitational delds, we use the statistical 


method of 


Chandrasekhar! (119431 ) . As we have shown above, within this method, the 


distribution function is never Gaussian for any long-range potential, obtained as multipole 
expansion of Newtonian one. At the same time, all previous approaches postulated (rather 
then derived) the distribution function in Gaussian form, which, in our opinion, does not 
redect the physical nature of long-range gravitational multipole interaction, which generates 
distribution functions with long tails. 


3.1. Numerical calculation of f{E). Dimensionless variables. 


Following 


Ghandrasekharl fjl943f) for the case of Holtzmark function, we introduce the 


dimensionless variables pE = x and P = E/a^^^. In these variables the integral ffTSj) renders 
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to 


fm = 


Jo 


X sinx exp 




3 / 4 ' 


Hm = 


ttP 


X sinx exp 


-'i) 


PJ 

3 / 4 ' 


... 


47r/3^a^ ’ 


dx. 


( 21 ) 

( 22 ) 


The physical meaning of the function H^fd) is that as it proportional to f3‘^f{(d), it is just the 
integrand in Eq. flT^ . being the effective one-dimensional distribution function of random 
helds. In other words, the normalization condition for H{f3) assumes one-dimensional form 


/ = 1 . 

Jo 

In this case, the average value of dimensionless random held fd reads 

poo 

j}= j j3H(j3)d^. 


(23) 


(24) 


Function H{ld) will be calculated numerically below. 


3.2. Asymptotics of distribution function f{(d) 

We begin with asymptotics of H{(d). At /9 —)■ 0 we make substitution x/(d = t to obtain 
from (12^ 


H{fd) = ^ / tsinfdt ^ ^ / t(pt- -{fdtf + 

TT Jo tt Jo V 6 J 


e 


(25) 


The hrst term of the expansion 

2/52 roo 


yields 


H{ld ^ 0) = 


TT 


Svr 


.2 -*3/4 . 8/12 p 3 _ 4g 

re dt = / re dz = —fd- 

Jo 


Stt 


(26) 
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At /3 ^ oo we expand the Eq. fl2^ over small parameter 1//1. We obtain in the hrst 
approximation 


2 r 

Hi(3 oo) = — lim / 


xsinx ( 1 


X 


3/4 


/ 33/4 


e-^^dx = 


2 

lim 




x'^/^sinxe ^^dx = 
r(ll/4) cos - 0.945972642/3-^/h 


TT/I’^A 5^-0 
2 


(27) 


To derive expression fl27|) we take into account that lim^^o x sin xe ^^dx = 0. 


Having asymptotics H{{3), we calculate those for /(/5) with the help of relation 

f3 ^0 


= 


(28) 


/3 ^ oo. 


Here a is given by Eq. flT8|) . The asymptotics fl28|) shows that f{(3) does not depend on 
jS at small (3 and decays at large (3. This shows that although normalization integral is 
convergent (we recollect that normalization condition for f{(3) looks like °/3^f(/3)d^ = l), 
already hrst moment does not exist. This can also be seen from large f3 asymptotics of 
//(/3) fl27D . The expression (128|) is a conhrmation of the fact that function /(/3) belongs to 
the class of heavy-tailed distributions. 


4. Distribution function of angular momenta 


Our aim is to calculate the distribution function of galaxies angular momenta Lj. For 
that we need a relation between the angular momentum L of a galaxy and its gravitational 
held Eq„ad(r) ([3]). The expression for components (a = x,y,z) of L has been derived 
in the form of perturbation serie s in small Lagrangian co ordinate q. The hrst order terms 


are dehned by Eq. (11) of Ref. fICatelan fc Theuns 


19961) . while second order ones by Eq. 
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(28) of the follow-up article fjCatelan fc Theuns 


1996ah . 


structure and can be written in the form 


Both expressions have identical 


r® = a, /3,7, ff = X, y, z, 


(29) 


where index i = 1,2 denotes the order of perturbation theory, ea/ 3 -y is Levi-Civita symbol, 
are components of quadrupole (tidal) held (jS]) and la^y are components of inertia tensor. 
Accordingly, functions fi(t) = a?{t)D{t) and / 2 (f) = E{t) (dot means time derivative) are 


known functions of time, calculated from t 


of perturbation theory flBouchet et al. 


le differential equations, derived in i - th order 


19921) . The explicit form of these equations read 


tlb{t) -h a{t)D(t) = 0, 
tlE{t) + a{t)E{t) = —a{t)D{tY, 


(30) 

(31) 


where 0 < t < 00 is dimensional physical time and to is some characteristic time, depending 
of the cosmological model considered. We will choose this quantity below. The dimensionless 
function a{t) (so-called scale factor) is determined from zero-order perturbative equation, 
which indeed is hrst Friedmann equation, depending again on chosen cosmological model. 
General form of this equation reads 

H‘2 

-—2 = Gijo ^ -l- ^ -1- G/jO ^ -|- Ga- (32) 

-^0 

Here El = a/a is Hubble parameter (a = da/dt), Ho is Hubble constant and Qi {i = 
R,M,k,A) are corresponding density parameters taken at present time, when a{t) = 1. 
Specihcally, Qr is radiation density, Qm is matter (dark plus baryonic) density, flk is 
co-called spatial curvature density and Ga is cosmological constant or vacuum density, 

Ga = where A is cosmological and Ho is Hubble constant. For our calculations of 

distribution function of angular momenta, we will choose hat ACDM model of the Universe, 
keeping in the equation fl32|) only Qm and Ga terms, Qm -I- Ga = 1- 
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Although the components of the gravitational tidal (shear) field E are different in the 
first and second orders of perturbation theory, for our purposes it is sufficient to consider 
them to be the same as they are simply the arguments of distribution function flTK]) . 
Denoting b(t) = fi{t) and omitting index i in the components of tidal field E, we can rewrite 
the relation explicitly: 


Ljb{t) 

Ly/b{t) 

Ljb{t) 


E I E I +E I -E I -E I -E I 

■‘-'yx-^-xz ~ -^yy-^-yz ' ■‘-^yz-^zz ■‘-^zx-^-xy ■‘-^zy-^yy ■‘-'zz-^-zyi 

E^ZX^XX E^ylyx H“ Exx^XZ Exylyz Exz^ZZ-) 

El -\- E T -\- E I —El —El —El 

■^xx-^xy ~ ^xy^yy \ ^xz^zy ^yx-^xx ^yy^yx ^yz^zx' 


(33) 


Taking into account the symmetry relations lab = ha and Eab = Eba and leaving only E^z-, 
we obtain from 033 p 


Lx Ly b{t^Ez:zIxzi L^ 0 . 

^Ll + Ll + Ll = b(t)E„^Il + 4^. = UE, (34) 

U = Ut) = + 

The expression (l34p constitutes linear relation between angular momentum and tidal field 
moduli both in linear [i =1) and nonlinear [i =2) regimes. 


As the relation between gravitational field modulus and angular moment fIMP is 
linear, the shape of distribution function of angular moments f{L) is similar to that of 
gravitational fields. The explicit transition from f{E) (|T^ to f{L) can be accomplished 
combining the expression fl34P and known relation from the theory of probability 


/(L) = f[E(L)] 


dE 

dZ 


(35) 


which yields 


f{L) = 


271^1 


P e 




Sin I p 


Lo{t) 


dp, 


(36) 










19 


where -ho(t) is given by expression flMj) . Passing to dimensionless variables 

L- 


(37) 


generates the following pair of functions similar to the case of gravitational helds distribution 

(in]), dnD 


/(A) = 


27r2A^a‘^Lo 


X sinx exp 


2 


X smx exp 


xx 3/4 


XX 3/4 

x) 


dx. 


dx = 


H{X) 


dvrA^a^Lo ’ 


(38) 


(39) 


In this case the effective ID distribution function H[X) is similar to that from the 
gravitational helds (12^ and we have the asymptotics (125]) (divided by Lq) for the 
distribution function of momenta. 

The effective ID distribution of helds or momenta is shown in Fig. [T] It is seen that 
this function has characteristic bell shape and is asymmetric. Asymptotics fl27|) shows that 
the integral, dehning the mean value of the galaxy orbital moment (I24|) . is divergent. To 
estimate the most probable value of the orbital moment, we calculate A^ax, corresponding 
to the maximum of distribution function H{X), see Fig. [T] for details. This situation is 
typical for so-called heavy-tailed distributions like Cauchy one f{x) ~ (a^ + x‘^)~^. Such 
distributions can be met very frequently in all branches of physics, dealing with random 
processe s and ranging from conden sed matter physics to chemical kinetics and econophysics. 


see, e.g. 


Kapur fc KesavanI ()l992l) . As Cauchy and many other heavy-tailed distributions 


do not admit a hrst moment (corresponding integral is divergent), the maximum of such 
probability density function is usually taken as a measure of its mean value. 


The analysis of this mean value will permit us to derive some useful relations, 
which earlier had been taken as empirical ones. To perform this analysis we adopt 












the simplest possible C DM cosmology in t 
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le first order of perturbation theory, where 


a{t) = D{t) = ( Doroshkevich 1970 ) so that Lq = r = t/to, I = 


The solution of the equation dH/dX = 0 reads 


XmnT = 0.602730263. 


(40) 


In dimensional units fl37|) we have from 


= 0.602730263 ■ (27r ■ 0.41807255)^/V/^L, 


GQ 


= 1.09228264 (—j LqGQ, 


= 0.7281884 


NX 


1 / 


4/3 




N\ 




4/3 


t 


grW 


(41) 

(42) 


where k is a constant of order unity. To derive the last equation fl42|) . we estimate (on 
the base of Eq. ([2])) both galaxy quadrupole moment Q and its mean inertia moment I 
as being proportional to mR^, where m is galaxy mass and R is its mean radius. If we 
estimate volume V as V = R^, we see that R cancels down in Eq. fl42|) so that we obtain 
Lmax ~ 

On the other hand, we can suppose that volume V is related to the galaxy cluster so 
that its value V = R\, wh ere Ra is a mean radius of the cluster, which is proportional to 


the autocorre 


19971 : 


Longair 


ation radius fjPeebles 


1973 


1980 


Peacock 


1999 


Lin et ah 


1996 


Tucker et al 


20081 ). Although n is still a constant for any particular cluster, it varies from 


cluster to cluster with increasing richness N. In this case we may rewrite N = M/m to 
obtain the different form of expression for Lmax 

4 


/ ! 

■‘-'max 




which does not contain p. It seems formally, that Eq. 


implies that 


(43) 
t, but 


the time evolution of galaxy radius R and mean cluster radius Ra, being very complex 
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astrophysical process flLongairll2008l) . can complicate real time dependence Lmaxit) a lot. 


This question needs to be studied additionally. 


The equation (H2D shows that mean orbital moment of a galaxy is proportional 
to where N is number of galaxies. The dependence Lmax{N) = {k 2 = 

(USD is shown in FigJ2J It is seen, that in our model of constant galaxies 
density n = const, the systems with larger number of galaxies have larger angular momenta. 
Below, comparing the linear and nonlinear regimes of fluctuation growth, we will show that 
the assumption of constant density may safely be used for qualitative analysis of angular 
momentum acquisition. 


Let us hnally pay attention that the dependence between mean momentum of galaxy 
cluster and its mass comes from parametrical dependence of L^ax on galaxy mass, volume 
and quadrupolar moment. This same dependence can be identically rewritten in different 
ways. Then, if we assume that certain parameters are constant, we obtain different 
dependences of Lmax on not only cluster mass M but on mass and galaxies density as well 
as on galaxies number N. For instance, in Eq. flT2]) . at i? = const L^ax scales as a square of 
galaxy mass m in contrast to relation fl43l) . This shows that observational verihcation of the 
dependences L^ax {N, m) would permit to make unambiguous conclusion about constancy 
of a particular stellar parameter. 


4.1. Time dependence of distribution function in ACDM model. 


The time evolution of distribution function fl39D can be obtained with respect to the 
dehnition of A flTTl) and subsequently Lq fiMl) . The dependence L^it) generates substitution 
A —)■ \/ fiir), where r = t/to so that we have from fl3^ 


/ x smx exp 

Jo 






2 

ttA 


dx, i = 1,2. 


(44) 
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To derive Eq. fHT|) we take into account that there is additional coefficient Lq in the 
denominator of fl55]) containing /i(r) so that there is no /i(r) before the integral in Eq. 
fHT|) . To obtain /i, 2 ('r) in ACDM model, we begin with the determination of a{t) from 
Friedmann equation (j32]), which reads 

^ (45) 

at y a 


The solution of the equation fl45|) has the form 


a{t) = asinh^^^(f/fo), a = 


1 — 
12a 


1/3 


) to — 


3Ho^/^A 


(46) 


Now, the function /i(r) = a^{T)D'{T) {D' = dD/dr), where D{t) can be found 
numerically from the equation fl30|) . Accordingly, in the nonlinear regime, the function 
/ 2 (r) = E'{t) should be found numerically from the equation (13T1) . 


We note here, that in Einstein - de Sitter m odel /i(r) = (2/3)r and / 2 (t) = (—4/7)r^/^ 


flDoroshkevich 


197nl: 


Catelan fc Thenns 


occurs at 


1996al ). Also, the maximum of the function 
XmaAr) ~ 0.602730263/,(r), (47) 


where Xmax ~ 0.602730263 is a maximum of time - independent function H{X) (|40|) . It 
is seen that functions / 2 (t), related to the second perturbative corrections, are negative. 
Substitution of such function into the exponent of integrand fl4T)) generates the imaginary 
part, which does not change the behaviour of H{X,t) qualitatively. That is why for the 
purpose of comparison of the linear and nonlinear regimes of fluctuation growth, everywhere 
we use the moduli of the functions / 2 (t) both in CDM and in ACDM models. 


The dependences H{X,t) fl44D for Einstein - de Sitter CDM model with above fiij) 
are shown in the Fig. [3l It is seen that at time growth the distribution function diminishes, 
while at zero time it goes to inhnity. As time (hgures near curves in Fig. [3]) grows, the 
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maximum Xmax{T) shifts towards larger A so that the whole distribution function ’’blurs” 
at large times. This is related to the fact that functions /^(t) enter the exponent in the 
integrand (jHj). The comparison of left and middle panels of Fig. [3] show that the behaviour 
of H{X, t) is qualitatively similar in linear and nonlinear regimes of flnctuation growth. 
This means that for qualitative analysis we may safely use the linear regime. To further 
demonstrate that, at the right panel of Fig. [31 we report the dependence Xmaxij) dlZD- If 
is seen that both in linear and nonlinear regime this function grows with time, although in 
ACDM model this growth is much faster so that to show both CDM and ACDM curves in 
one panel, we use the logarithmic scale. 


The result of numerical calculations in ACDM model are reported in Fig. 01 It is seen 
that qualitative behavior of FT (A, r) is similar to that for CDM model. Here, however, 
we can trace the variation of H{X,t) with parameter a. It is seen that as parameter a 
increases, the distribution function decreases - for the similar times the distribution function 
is smaller for larger a. Similar to CDM model, the maximum of distribution function 
shifts towards larger A at time growth. The above regularities are the same for linear and 
nonlinear regimes, while the values of distribution functions in ACDM model for the same a 
are smaller in nonlinear regime. This is related to the fact that function / 2 (t) grows much 
faster then fi (t) in ACD M model. Also, as has been noticed in many references (see, e.g. 


Catelan fc Theunsl fjl996l) and references therein), that first order perturbation result (linear 


regime) corresponds to so-called Zeldovich approximation which is approximately valid also 
for nonlinear sitnation. This shows once more that for qualitative discussion of the time 
dependence of the distribution function H{X, r) we can safely use the linear approximation, 
based on the function /i(r) only. It should be noted that parameter to in Eq. fi46|) is of order 
of the Hubble time I/Hq. The comparison of the curves for different t = t/to in the Fig. 0] 
shows that the relaxation time is very long. This is in agreement with the fact that clusters 
are known to be dynamically young objects, i.e. with crossing timescales which are non 





24 


negligible with respect to the age of the Universe at the time of their formation. The fact 
that the relaxation time in Fig. |l]is very long snggests that the dependence of the effective 
ID distribntion fnnction H{X,t) (and average angnlar momentnm as the result) on the 
redshift should be weak. Also, the comparison of Figs [3] and |4] shows that the distribution 
of spin parameters of galaxy clusters depends on the cosmological parameters and generally 
speaking is weak. 


5. Discussion. Relation to observational resnlts 


Our theoretical results about the distribution function of momenta can shed some 
light on the problem of galaxies orientation in clusters. Our main message is that although 
the gravitational interaction between galaxies is of long-range character, the observations 
(which we will discuss below) may evidence that there is additional short-range intergalaxy 
interaction with characteristic radius Vc- This means that if the distance r between two 
galaxies is smaller then r^, they are correlated and have their orbital moments aligned. This 
is the case for the dense (rich) galaxy clusters, which, by this virtue, have high degree of 
orbital moments alignment. In the opposite situation of poor clusters, where the intergalaxy 
distance r > Tc, the long-range multipole interaction of alternating sign dominates and the 
alignment of the orbital moments is absent. We speculate that this situation resembles 
that in diluted magnetic systems, where the presence of short-range exchange interaction 
between magnetic spins promot e long-range magnet i c order, which is character i zed b y 


macroscopic spin alignment, see 
references therein for details. 


Stephanovichl (119971) ; ISemenov fc Stephanovichl (120031) and 


The aforementioned statistical method permits to account for this situation if we add 
the (empirical) short-range interaction term to the initial potential ([3]). In this case, the 
distribution function of random helds would depend on the above average value of angular 
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mome ntum Lmax = Lav as a parameter fsee IStephanovichl fll997l) : ISemenov fc Stephanovich 


(120031) ) so that self-consistent equation for Lav of tho form 


/ L(E)f(E,L,,)d^E 


(48) 


can be derived. Here f{E,Lav) is the distribution function of random gravitational helds, 
which substitutes the expression (1T8|) in the case of inclusion of the possible short-range 
interaction term. In such case, for hnite r^, the distribution function decays at if —?■ cxd 
faster then (IlSp so that the integral (I48p converges. As now the total interaction potential 
includes both luminous and dark matter, the equation (HHP permits to address the question 
about alignment of sub-dominant galaxies, when most of cluster angular momentum 
is in the smooth dark matter halo component. Observationally it is relatively easy to 
analyse the orientation of angular momenta in the luminous matter i.e. in real galaxies 
and their clusters. With dark matter (sub) halos this is not so easy. One should not. 


however, forget, that there are relation between the 


matter. The articles fjPaz et ah 


2008 


Bett et ah 


properties o 


2010 


Kim et ah 


luminous and dark 


2011 


Varela et ah 


201211 


show clear observational evidence of re l ation between the dark matter halos and galaxies 
orientation, see however 


Tenneti et ah 


(1201511 for opposite opinion. The assumption that 


angular momentum of the luminous matter traces that of the associated dark matter 
(sub)-structures (e.g. angular momentum totally aligned, one galaxy per sub-halo) allows 
us to conclude that angular momentum alignment of galaxies give us information about 
similar alignment in dark matter. In other words, the properties of angular momentum 
of luminous matter (like real galaxies) give us information about those of dark matter 
(sub)-structures. 


For luminous matter, it is possible to consider the relation between 


mornentum and the ma ss of a structure as based on the observational data (IGodlowski et ah 


2005 


Godlowski 


he angular 


2011 aJ). It is possible to investigate how this image varies with the mass of 
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galaxy clusters, beginning with the simplest ones, i.e. galaxies pairs. These investigations 


'Karachentsev & Mineva 

1984a 1 

): Mineva 

1987) 


Heloul fjl984l) examining a sample of 31 


galaxy pairs, found that an ’’anti alignment” of these galaxies spins occurs. 


Parnovskv et ah 


(119971 ) recognized a weak alignment in physical pairs of 


spiral galaxies was also discerned by 


g alaxie s. Alignment in pairs of 


in galaxies pairs has independently been conhrmed by 


Pestana. fc Cabreral (l2004f). Intrinsic sp in alignment 


Hevmans et al. 


( 2004 ) within their 


research on weak gravitational lensing, where it was necessary to estimate and remove the 
effects related to alignment of galaxy orientations. Also the analysis of positions of the Milky 
Way’s companions shows their non random distribution (they are located perpendicularly 
to the Milky Way’s disc), which can be regarded as their orbital al ignment. Galaxies withi n 


compact groups rotate on prolate orbits about the group’s c entre (ITovmassian et al. 


which contribute s to th e system’s to 
' Wang et al. 


Sales fc LambasI (1200 


al angular m omentum. (lYang et al 


(12009 


20011 ) . 


20061) found, while 


20101) conhrmed it, that the companions of 


central red gah 

axies are ah 

med along their large axes. The similar result has been obtained 

by 

Ibata et al. 

(2013, 

2014 

) in the two papers about the ordering of satellites orbits around 

M31 Theirs latest article ( 

Ibata et al. 

2015 

) also corroborates this result. Thus it can be 


maintained that structures like galaxies and their companions, pairs and compact groups 
of galaxies have a non zero net angular momentum related mostly to their orbital motion. 
One should not forget, however, that situation is more complicated in more massive 


structures. As there is no evidence for rotation of t 


re groups and clu s ters of galax i es (se e, 


e.g.. 

Regos & Geller 

(1989) 

Diaferio & Geller 

(1997 

); 

Diaferio 

(1999) 

Rines et al. 

(2003) 

Hwang & Lee 

(2007’ 

), it is clear that angular momentum of such structures is related 


primarily to the alignment of constituents spins. 


We should note here that there are many (seemingly contradictory) observational 
results regarding the alignment (or misalignment) of galaxies angular momenta in the 
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literature beginning with the paper of iThompsonI fjl976l) . who found an alignment in 
the galaxy orientations in the Virgo and A2197 clusters. Th e evidence for alignment 


MacGillivrav& Dodd 

(1985a). Non random ealaxv orientation 

was found in very rich galaxy 

1 

clusters (Diorerovski! 1 

983; 

Godlowski et al. 199 

8; Wti 

1998 

: Baier. Godlowski & MacGillivrav 

2003: Kitzbichler& Saurei 

2003: 

Arval! 

2007: 

Godlowski 

2012 

). The alignment of galaxy 

planes was also found in A1689 

(HnnET et al. 

Mio[ 

dune & Ebeline 2012). This result is 


important as A1689 is the most distant cluster where the alignment has been found till now 
Th ere are also co n tradic tory results. For instance, 


(12Q03I); 


Hofman at ah 


Bukhari! fll988l) : 


Bukhari fc Cram 


(11989!) studied orientation of galaxies wit hin clusters and did 


not hnd any alignment. The same result has been obtained by 


Arval fc Saurer! ( 20051) 


in theinvestigations of three Abell clusters of ric h ness class zero. D u ring studie s 


of the isolatedAbell groups (iFlin fc Olowin 


Niederste-Ostholt et ah 


1991: 

Trevese et al. 

1992; 

Kim 

2001 


20101) . only a rudimentary alignment was found and related only 


to the brightest cluster members. The alignment has also not bee n found during analysis 


of Tu 

Iv’s eronps of ea 

axies beloneine 

to the 

1999; 

Godlowski et al. 

2005; 

Godlowski 

2011h) 


o the Local Supercluster (iGodlowski fc Ostrowski 


2011bf) . Summarizing above observations of angular 


momenta misalignment, we can conclude that such misalignment had been obtained for less 
massive structures like group and poor galaxy clusters. 

To check the hy pothesis, that the alignm ent of galaxies angular momenta increases with 


the cluster richness. 


Godlowski et al 


(120101) examined orientation of galaxies in clusters 


both qualitatively and quantitatively. The analysis of the spatial orientations of galaxies 
in the 247 optically selected rich Abell clusters, h a ving i n the considered area at least 100 


members has been performec 


from the PF catalog (jPanko 




Godlowski et al. 


(120101 ). The structures have been taken 


20061) . The statistical analysis, based on linear regression. 
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permitted to conclude that cluster angular momenta increase with their nu merousness. Note 
howev er, that relatively small statistical sample of 247 clusters, analyzed bv iGodlowski et al. 


11201011 ■ does not give a possibility to discriminate bet ween linear dependence te sted by 
Godlowski et al.l ( 2010 1. the dependence Lmax ~ f Gatelan fc Theuns ( 1996 11 and our 


result (l43l) that mean angular momentum is proportional to (stemming from 

However, such detailed analysis would be possible if larger statistical sample of galaxy 
clusters is available. 


The above results show clearly that galaxy clusters have a non zero net angular 
momentum related mostly to their orbital motion. For more massive structures there is a 


lack of alignment of the orientation of galaxies for group and poor galaxy clusters, w 


rile 


ther e is evidenc e for al ignment for the rich clusters of galaxies (IGodlowski et al. 


also 


(1299511 . see 


Godlowski! (j2911a| l for later improved analysis). It has been suggested that degree of 


alignment increase with clusters richness, and as a result the cluster angular momentum 
increases with its numerousness. Here we emphasize once more that the equation (148|) . 
which depends parametrically on the clusters parameters (like mass density p and cluster 
density n) will give zero or nonzero solutions for Lav (corresponding to alignment or 
misalignment of galaxies angular momenta) depending on the parameters values and cluster 
reachness in particular. 


Now we present an alternative explanation of alignment and misalignment of galaxies 
angular momenta in rich and poor clusters respectively. Namely, as rich galaxy clusters 
are less dynamically evolved objects, the galaxies constituting them, haven’t had time to 
reach the pericenter of their orbit in the cluster and retain the alignment imprinted by 
the larger scale environment in which the cluster is embedded. In other words, large scale 
hlaments (cosmic web) feed (reach) galaxy clusters along privileged directions, resulting 
in clusters being rather more prolate in shape than spherical. We begin with the paper 
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of iGodlowski et al.l (119981) . who had shown that the galaxies orientation distribution in 
the Abell 754 double cluster is nonrandom with galaxy planes being perpendicular to 


the main cluster plane. The a bove nonrandom orientation distri 


Dution had also been 


found in the Abell 754 cluster (iBaier. Godlowski fc MacGillivrav 


20031) . but the direction 


(relatively to the main cluster plane) of the observed galaxy ordering is perpendicular to 


that for Abell 754. T 


re interpretation of above orientation difference has been presented by 


di Fazio fc FlinI (119881) in terms of two different types of galaxy clusters: oblate and prolate. 


One more interpretat ion can be done on the base of tidal interaction scenario. Namely, it 
has been observed by 


Paz et al. 


(120081) that in large scale structure the direction of angular 


momenta (relatively to its main plane) of co nstituting objec t s dep e nds on the s 


mass. The same resu' 


he n e west analys i s (IGodis et al. 


( 2000 ): 


(120071) : 


Lee fc Pen 


Paz et al 


t has been obtain ed by 


(12000 


20081) : 


Truiillio et al. 


(120061) : 


Varela et al 


2012) (based on t 


le earlier studies of 


Bailin &: Steinmetz ( 2005 ): 


Zhang et al. 


Aragon-Galvol ( 20071) : 


ructure 


( 120121 ). 


Sugerman et al 


Hahn et al 


(120091) ) on the dark matter halos angular momentum 
orientation also conhrm the above result. Th e studi e s of galaxies angula r momenta ordering 


in large scale had been fulfilled by 


Paz et al 


(I2OO8I) : 


Zhang et al. 


(I2OI3I ). who use the data 


from Sloan Digital Sky Survey catalog. They found that galaxies angular momenta align 
perpendicularly to the large scale structures planes. Latter effect has not been observed for 


the structures w ith relatively small masses. These result s agree well with 


:he simulations of 


Paz et al 


(120081) based on tidal interaction mechanism. 


Jones et al. 


( 2010 ) have found that 


the spins of spiral galaxies in the cosmic web have tendency to align along the filaments 
axes, which has been interpreted as the ’’fossil” evidence of the effects of long-range tidal 
interactions. 

The other interesting problem is possible time evoluti on of galaxy clusters alignment. 


Assuming Einstein - De Sitter model (jPoroshkevich 


19701 ). on the Fig. |3]we report the 


time evolution of the function H{X). It is interesting that at time growth the distribution 
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function goes to zero. This means that older structures (clusters) should have more 
scattered observed values of angular momenta than younger ones. At the same time, 
equation shows that mean angular momentum of the clusters should increase in time. 


Latter result is obtained in Einstein - De Sitter model but is still valid for any similar 
cosmological models. The predictions will be available to verify when we get better data 
concerning alignment in galaxy clusters. Note however that even now we have observational 
results sugges ting the pos s ible ev olution of alignment with redshift. There are for example 


the results of 


Song fc Leel (120121 ) who found that the alignment prohle of cluster galaxies 


drops faster at higher redshifts. 


6. Conclusions 

The main physical result of the present paper is the calculation of the distribution 
function of the gravitational fields of astronomical objects like galaxies ensembles and 
smooth halos based solely on the tidal interaction between constituting elements. We 
show that for tidal (quadrupolar) interaction the distribution function cannot be Gaussian, 
its explicit form being presented by the equation (ITSll . We emphasize here that derived 
distribution function of gravitational helds does not depend on the specific Eulerian or 
Lagrangian description of Newtonian matter and thus can be used to calculate virtually any 
observable characteristic of stellar ensemble. As an example, we use the above distribution 
function to calculate the distribution of angular momenta. From the astronomical 
interpretation point of view, it is important that for particular cluster with richness n we 
expect not a single value of angular momentum but the range of allowed values described 
by the probability function. As the distribution function ([15]) slowly decays at infinity, its 
first moment does not exist because the corresponding integral is divergent. To calculate 
the mean value of angular momentum in this situation, we assume that the maximal value 
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of distribution function Lmax gives the desired quantity. We note that such procedure is 


usual while dealing with so-called heavy-tailed distributions, see, e.g. 


Kapur fc Kesavan 


fll992l) and references therein. 


As astronomical objects (for instance galaxies) masses, radii and number enter the 
equation fllSp as parameters, we were able to show that mean value of angular momentum 
Lmax for partic ular galaxy clust er of mass M is proportional to thus corroborating 


well-known (see ISchaeferl (120091) and references therein) empirical result. Our other result 

LifjlClX N^/3 ( 

or eq. (143|) ) has also its astronomical interpretation that larger (richer) 

clusters of galaxies have higher angular momentum, see also the discussion below. The 
observational discrimination between the dependences and would be possible 

when larger statistical sample of galaxy clusters will be available. The parametric time 
dependence of Lmax via functions a{t) and D{t) in Einstein - de Sitter model permits us to 
trace its time evolution. This shows that our approach to derivation of distribution function 
of galaxies angular momenta gives physically reasonable answers. We have also analyzed 
the time evolution of the distribution function. It is reported in Fig. |3]and shows that the 
distribution function is flattening with time. 

The relation between angular momentum anc 


extensively analysed theoretically 


1982 


Sistero 


1983 


Brosche 


1986; 


(iMuradvan 


1975 


Mackrossan 


19871: 


mass o: 

the st 

ructures has been 

Wesson 

1979 


981 


1963; 

Carrasco 

Paz et ah 

200^ 

1), see 

Schaefer 

(2009) 


for review. This relation has usually been presented empirically in the form Lmax 
From the point of view of modern scenarios of galaxy and their structures formation. 


increasing of angular momentum with the cluster richness could be ex 


torque scenario in the h ierarchic clustering mo del (iHeavens fc Peacock 


2007 


Noh&Lee 


2006alJbl ) and in Li model ([L 


jlained only in tidal 


1998 


Godlowski et al 


note however that the value of the Universe rotation, required by 


1988 


Hwang fc Lee 


2005h . One should 


Lil (119981) . is too large as 
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compared to the anisotropy found in cosmic microwave background radiation (CMBR). 


The increase of the galaxies angular momentum with mass of the s tructure was foun d 
observationally during analysis of t 


2005 


Arval 


2007 


Godlowski et ah 


re alignment of galaxies i n clus ters (iGodlowski et ah 


2010 


Godlowski 


2011a 


agreed that grou 


3S anc 


Diaferio fc Geller 


1997 


clusters of galaxies do not ro t ate flRegos fc Geller 


Diaferio 


1999 


Rines et ah 


2003 


20121). Since it is com monly 


Hwang fc Lee 


1989 


20071) . any possible 


nonzero angular momentum of such structures should arise only from possible alignment of 
galaxy spins and stronger alignment mean greater angular momentum of such structures. 
Generally there is no evidence for a non zero angular momentum of galaxies or group 
and poor galaxy clusters, while we observe such alignment for rich clusters of galaxies 
and superclusters. We speculate that this phenomenon may occur due to some additional 
short-range interaction between galaxies such that in rich clusters the galaxies are correlated 
as they fall in the range of this interaction and hence have their angular moments aligned. 
In such situations there should be some critical richness ricr, related to the interaction 
range Tc such that only clusters of richness n > Ucr would have their spins aligned. We note 


here that sue 


Stephanovich 


1 phys ic al picture is common for disorde red magnets and ferroelectrics, see 


(119971 ): ISemenov fc Stephanovichl (120031) and references therein. We postpone 


the quantitative investigation of this interesting question for future publications. The above 
scenario can we applied to the problem of the possible merger of galaxies into larger objects. 
The corresponding results will also be published elsewhere. 


The problem of galaxies merger in a cluster is related to that of the role of much more 
massive central galaxy. This problem, in turn, is due to the fact of (generally speaking 
random) interaction between cluster members and dynamic evolution of the nearby (to 
specihe cluster member) structures. The problem of dynamic evolution can be studied 
by the combined examination of mutual orientation of galaxies in clusters and Binggeli 
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effect ( 

Bineeeli 

1982 

. c 

Itruble & PeebleJ 

1985; 

Ghambers et al. 

2000 ; 

Hashimoto et al. 

2008 

Godlowski & Flin 

2011 

1). More specihcally, this can be done by two methods. The hrst is 


Binggeli effect studies, i.e. the investigation of relation between positions of great axes in 
groups or clusters of galaxies and directions towards their neighbors. The second one is 
the studies of mutual orientation of the brightest galaxy (and other bright galaxies) in a 
structure relatively to the position of cluster great axes or even examination of structure 
ellipticity redshift dependence, especially in the enlarged observational samples. The 


analysis of the differences betwee n position angles o 


to the Local Supercluster shows flGodlowski fc Flin 


the T ullys groups of galaxies belonging 


2 OIOI) that there exists the alignment 


of the line joining two brightest galaxies with both the position angle of the parent group 
and the direction toward Virgo cluster center. This analysis reveals the following picture of 
the structure formation. Two brightest (most massive) galaxies were formed hrstly. They 
originated in the hlamentary structure directed toward the center of the protocluster. This 
is the place where the Virgo cluster center is located now. Due to gravitational clustering, 
the groups were formed in such a manner that galaxies follow the line determined by the 
two brightest (most massive) objects. Therefore, the alignment of the structure position 
angle and line joining two brightest galaxies is observed. The other groups were formed on 
the same or nearby hlament. This shows the particular role played by the more massive 
(brighter) galaxies. In our future investigations, we will analyze the particular role played 
by the central (most massive) galaxies quantitatively. 

One should note that the theoretical and observational analysis of galaxies alignment 

point of view of weak gravitational leasing 


(Troxel & Ishak 

2014 

Joachimi et al. 

2015) 


2015h . There are mutual influence of the orientation 


of galaxies and weak gravitational lensing. It should be pointed out as well that the 
examination of galaxies orientations is also meaningful due to one of the ou tcomes of the 


activity of gravitational lensing effect flHeavens et al 


2000 


Schneider 


2005 ) which is the 
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alignment of the galaxies imag es. Such alignment, in s everal Mpc scale is also expected in 


case of cosmic shear existence. 


Crittenden et al 


(1200 ll) proved that at least in the scenario 


of tidal interactions, the effects of alignment can be distinguished from the effects of 
weak gravitational lensing. Taking both of these effects into consideration (in appropriate 
proportions) is of a paramount importance for mapping the mass distribution with weak 
lensing techniques, and vice versa: weak lensing induced shape deformations are important 
for studies of the intrinsic orientation of galaxies within structures. Therefore weak-lensing 
studies will allow investigating mass distribution in clusters which is important for studies 
of dark matter in them. 

Let us hnally summarize the simplifying assumptions made in the above calculation 


of the distribution function. We assume t 

hat autocorre 

ation radius is constant and is the 

same for each cluster ( 

Peebles 

1973, 

1980; 

Peacock 

1999; 

Lin et al. 

1996; 

Tucker et al. 

1997 


Longair 


200811 . Also we treated cluster density as a parameter rather then a function of 


interobject (intergalaxy) distance. We note here that as time (via the functions 
calculated in CDM and ACDM models in linear and nonlinear fluctuation growth regimes) 
enters the distribution function as a parameter, our a pproach is valid for any cosmological 


model - from conventional Einstein - De Sitter CDM flDoroshkevich 


1970 1 to ACDM. Most 


important simplihcation is to assume that all galaxy has equal masses. In the future studies 
we plan to consider more realistic situation, introducing n = 'n.(r), the real (i.e. extracted 
from observational data) galaxy mass distribution function and the short-range term in 
the potential of interaction between astronomical objects (say galaxies and their clusters). 
Such generalisations will require numerical calculations of the distribution function of 
gravitational helds and its mean value. 
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Fig. 1.— The effective ID distribution function (13^ of dimensionless orbital moments H{X). 
The distribution function of dimensionless gravitational fields H{f3) fj2^ has the same shape. 
Dashed line shows the argument Xmax iHOl) . corresponding to maximum of H{X). Inset depicts 
real 3D distribution function 47ra^Lo/(A) = H{X)/X^. 
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Fig. 2.— Dependence of average galaxy angular moment Lmax = on number of 

galaxies in a stellar system. K 2 = . 
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Fig. 3.— Time evolution of effective ID distribution function H{X,t) in CDM model. Left 
panel - linear regime Xmax = 0.6027/i(r)=0.4018r, middle panel - nonlinear regime Xmax = 
0.6027 |/ 2 (r)| = 0.3444x^/^ (as / 2 (t) < 0, everywhere here we take its modulus). Figures 
near curves correspond to dimensionless time r = t/to. Right panel reports the dependence 
Xmaxij) for CDM and ACDM models. Solid lines correspond to linear regime, dashed - to 
nonlinear. Mind the logarithmic scale on the vertical axis, needed to place the curves of 
CDM and ACDM models having very different growth rates. For ACDM model we choose 

1/3 














Function H(?i, x), nonlinear regime Function H(A-, t), linear regime 


-48 - 




Fig. 4.— Time evolution of effective ID distribution function H{X,t) in ACDM model for 
different a. Upper panels - linear regime, lower panels - nonlinear regime. The notations are 
the same as in Fig. [3l 









